SPECTRAL CURVES AND PARAMETERIZATIONS OF 
A DISCRETE INTEGRABLE 3-DIMENSIONAL MODEL 



S. PAKULIAKj AND S. SERGEEVft 

Abstract. We consider a discrete classical integrable model on the 3-dimensional cubic 
lattice. The solutions of this model can be used to parameterize the Boltzmann weights of 
f~>) ' the different 3-dimensional spin models. We have found the general solution of this model 

constructed in terms of the theta-functions defined on an arbitrary compact algebraic curve. 
The imposing of the periodic boundary conditions fixes the algebraic curve. We have shown 
that in this case the curve coincides with the spectral one of the auxiliary linear problem. 
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In the case when the curve is a rational one, the soliton solutions have been constructed. 



GO 

1. Introduction 

The paper is devoted to a description of the periodic and soliton solutions of some generic 
classical 3-dimensional discrete integrable model. We will see that this description can be 
presented in the completely analogous way as the description of the finite-gap solutions of 
the hierarchy of the continuous integrable equations, although the way we get the system of 
equation is unusual and motivated by the approach developed in Q for discrete integrable 
spin models. First, we construct the discrete equations of motion from an equivalence of the 
linear systems (which replace zero-curvature condition for Lax operators in usual formula- 
tion) and after that prove the integrability by the counting the independent number of the 
"£h \ integrals of motion. 

j> | Several types of the boundary conditions may be considered on the cubic lattice (the open 

boundary, periodical boundary conditions in chosen directions, or the completely periodical 
boundary conditions). This choice leads to several dynamical interpretation of the model: 
as a Cauchy problem, as a Backhand transformation, or as an analogue of the standing 
vibrations on the discrete 3d torus. 

Starting from the discrete equations of motion we introduce a change of variables through 
triple of the Legendre variables, which transforms the equations of motion into the three- 
linear form. These three-linear equations appear to be a generalization of the famous Hirota 
bi-linear discrete equation. Then we observe that these three-linear equations can be formally 
solved with the help of Fay's identity for a theta-function on an arbitrary algebraic curve. 

Some facts observed in this paper are the demonstration of the general statements proved 
more than two decades ago in [Q] - any discrete integrable system can be solved using the 
algebraic-geometry methods. In this paper we develop further an alternative approach to 
3-dimensional discrete integrable systems |]J which does not use the notion of Lax operators 
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FIGURE 1 . Example of cubic lattice of the size 3x2x4. 

and can be applied equivalently to quantum (spin) or classical integrable systems associated 
with several 3-dimensional lattices. Instead of the Lax operators we use the notion of the 
linear system defined on auxiliary planes. In this paper we consider the cubic lattice, al- 
though these methods may be applied to arbitrary 3-dimensional lattices formed by a set of 
planes. 

2. Classical discrete integrable system on the cubic lattice 

Let the vertices of the cubic lattice spanned by orthogonal basis ei, e 2 and e3 be marked 
by the vector 

(1) n = mei + n 2 e 2 + n s e 3 . 

The cubic lattice is formed by three sets of parallel planes or, equivalently, by three sets 
of parallel lines (see Fig. [!]). Let us associate to each edge of the cubic lattice of a given 
type a, (a = 1,2,3 corresponds to three orthogonal directions of the lattice) the dynamical 
variables as shown in Fig. |2|. Namely, the pairs of dynamical variables u a>n , w a ^ n , a = 1, 2, 3, 
are associated with the edges incoming to the oriented vertex n, while u a ^ +ea , w a ^ n+ea are 
associated to the outgoing edges. Fig. |2| shows also two auxiliary planes intersecting the 
incoming and outgoing edges. Each of these planes cuts seven of eight octants around the 
vertex of the cubic lattice numbered by the vector n. 

Let us stress that our considerations at the moment are local. Our goal is to obtain the 
relations between dynamical variables surrounding given vertex. After this we extend these 
relations to the whole lattice and so obtain the discrete dynamical system. 

Each of the triangles on the auxiliary planes is formed by three lines obtained by the 
intersection of this plane with planes forming the vertex n of the cubic lattice. In this way 
the vertices of the auxiliary triangles are associated to the edges of the cubic lattice and so 
to a pair of the dynamical variables. Let us consider two linear problems attached to the 
auxiliary triangles via the following rules. First, we introduce the linear variables $ a , $6, $ c 
and $d around the vertex on the auxiliary plane according to Fig. |3|. The linear problems 
for each vertex on the auxiliary planes will be always of the form 

(2) = $ a - $ b • u + $ c • w + $ d • K u w . 
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FIGURE 3. The linear problem for the vertex with associated dynamical variables and 
additional parameter k. 



Here k G C is an additional parameter associated with the line of the cubic lattice. Note 
that the coefficients of the linear form @ are fixed by the orientation of the lines which form 
the vertex on the auxiliary plane. 

The enumeration of the linear variables on the auxiliary triangles is shown in Fig. It is 
clear that these linear variables are associated with internal part of the cubes in the cubic 
lattice. According to (Q) we can write the system of the linear equations 

^n+e2+e3^1,n ~l~ ^n^'l.n 4" ^n+es ^ l:n^l,n'^'l,n j 
= $ n - $ n +e 3 ^2,n + $n+e^2,n + $ n+ e 1 +e 3 «2 ;I1 «2, n W2, n , 
= $n+e 2 - $ n^3,n + &n+e 1 +e 2 w 3,n + $n+ei ^n^n^n 
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FIGURE 4. Parameterizations of the classical linear variables <& n in cubic geometry. 



for left triangle shown in Fig. |] and the linear system 



(4) 









^n+ei+62 ^ > n+e:L+e2+e3'U'l ) n+ei ,n+ei " 

^ ) n+ei+e3^1,n^l,n+ei^'l,n+ei i 
< ^ ) n+e2 < ^ > n+e2+es^2,n+e2 < ^ > n+ei+e2^2,n+e2~l _ 
+ < ^ > n+ei+e 2 +e 3 K 2,n'W2,n+e2' u; 2,n+e2 j 



n+e 2 +e 3 



n+e3^3,n+e3 



+ ( =&n+ei+e 2 +e 3 'U ; 3 jn+ e3 + 



^n+e2+e3 ,n+e3 

W 3 

,n+e3 



for the right one. 

Supposing that dynamical variables do not vanish identically on any edge of the lattice we 
may exclude from the systems © and (H) the linear variables $ n and $ n+ei+e2+e3 . Let us 
require that the systems (|3[) and (f|) (each contains now only two linear relations for the rest 
6 linear variables $ n +ei, $n+ei+e 2 5 etc.) are equivalent. This means that the matrix elements 
of these linear systems which are rational functions of all dynamical variables around given 
vertex coincide identically. Moreover, let us require also that each type of the parameters 
K ajn is conserved along the corresponding direction (which is already taken into account in 

(§) 



(5) 



^a,n+e a i 



a 



1,2,3. 



As the result, we obtain the recurrent relations for dynamical variables « Q ^ n , w a>ll1 u a)I1+eo 
and w n „i P „ : 



(6) 



u l,n+ei 



K2,nMl,n «2,n »2,n 



^l,n«l,n^2,n + «3,n «2,n «"3,n + «4,n «3,n Wl.n ™3,i 
W lja W 2 , D + U 3yn W 2 ,n + K3,n«3,nW3,n 



U! 3 ,x 



(7) 



U-2 
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,n+e 2 



^2.n- 



e 2 



«2,n«3,n + «2,n U>l,n + «l,n «l,n ™l,n 

Wl,n W 2 , n W3,n 

Wl,nW 2 , n + «3,nW2,n + «3,n W3,n ™3,n 



(8) 



«3. 



n+e 3 



e3 



f2,n«3,n + «2,nWl,n + Kl.n «l,n Wl.n 

«2,n«2,n^2,n W 3 , n 

«l,nWl,nW2,n + «3,n «2,n »3,n + «l,n «3,n «l,n ™3,n 



where due to (Bl) we have 



(9) 



K l,n — K l:n2,".3 ) K 2,n — ^2: ni 



ri3 j 



^3,n ^3:ni, 



ri2 • 



One may ask why the recurrent relations ©— © are the dynamical system. Let us introduce 
the following Poisson bracket for the edges incoming to the vertex with some fixed n: 



(10) 



{u c 



.} 



u n _ n w. 



a,n n 



and any other bracket for incoming edges vanishes. The transformation 
ical one, i.e. from (p~0|) it follows that 



is the canon- 



(11) 



and any other bracket for outgoing edges vanishes. In general, the auxiliary triangles shown 
in Fig. (Q) are the fragments of two adjacent space-like surfaces, and the set of 3d vertices 
between these surfaces provides by means of the canonical transformation of the set 

of dynamical pairs from incoming surface to the outgoing one. We call this transformation 
the local equation of motion in the direction perpendicular to the chosen space-like surface. 
Besides the local equations of motion a complete formulation of a dynamical system 

needs the specification of boundary conditions. 

The type of the boundary conditions as well as the global characteristics of the model de- 
pend on the choice of the space-like surface. The auxiliary lattices described are appropriate 
for the Cauchy problem for the cubic lattice. In Section |6| we will choose ni as the discrete 
time, while ri2 and n 3 will be the discrete space coordinates. At this choice the system (|6]) 
will describe the "time" evolution of the variables Ui >n , w^ n and the systems (0)-(H) - the 
space distribution of the "auxiliary" for this evolution variables w 2 , n , w 2,m u 3<n and u> 3 , n . 
All the dynamical quantities, such as integrals of motion, spectral curves etc., should be 
calculated in this case for the system (Ft). This will be subject of the Section H. 
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3. Legendre transform 

It follows from ©-© 

Wl,„W2,„ = Wi in+ei U>2,n+e 2 , U 2 , n U 3 ,n = M2,n+e 2 «3,n+e 3 , 

(12) 

u l,n M l,n+ei 



W 3 ,n ^3,n+e 3 

Relations (|T2|) provide in general the possibility of the following change of the variables: 

_ (0) r 2,n _ (0) 7"3,n+e 2 

M l,n — u l:n 2 ,n 3 > w l,n ~ w l:n 2 ,n 3 j 

r 2,n+e 3 73 )n 

. , (0) Tl,n (0) r 3,n 

(13) M2,n = *4ni,n 3 - > W 2,n = Wfcni.na ~ > 

'l,n+e 3 '3,n+ei 

(0) r l,n+e 2 (0) r 2,n 

«3,n = «3:m,n 2 - > W 3,n = ■ 

r l,n r 2,n+ei 

Eqs. ©-(U) with the substitution fll3|) become 



(14) 



^a,n'^a,n+e ( 3+e 7 '?73,n'7~7,n '?~a,n'77?,n+e 7 ^7,n+en ~r 



where (a,/?, 7) is any cyclic permutation of the indices (1,2,3), and the coefficients in (|T4| ) 
are 

(0) (0) (0) 

K 3:ni,n 2 W 3:ni,n 2 U l:n 2 ,n 3 U 3:ni,n 2 

S l,n : : (0) ' ri ' n ' 1 ~~ (0) (0) ' 

W 2:ni,n 3 W l: n 2 ,n 3 M 2: m,n 3 

(0) (0) (0) 

^3: ni,n 2 K 2: 711,713 M 2: 711,713 W 2: ni,n 3 



/i r\ d: 711,712 

(!5) s 2 , n = —7-: , r 2 , n 



(0) ' z > n (0) (0) 

1:712,713 I:n2,n3 3:ni,n2 1: 722,713 uy 3:ni ) ri2 



(0) (0) (0) 

"2:711,713 ^1:712,713^3:711,712 



S3 ' n (0) ' r 3,n ■ (0) (0) 

K l:n 2 ,n 3 U l:n 2 ,n 3 W 2:m,n 3 U 3:ni,n 2 

This change of variables, T a xi instead of u a>n and iv^n, has the following interpretation. 
Equations are the kind of the Hamiltonian equations of motion for the classical 

discrete system. Substitution ([13]) therefore is the Legendre transformation. And finally 
equations ([TJ|) are the Lagrangian equations of motion. We will call r a n as the Legendre 
variables, and the coefficients Ua:n ,n^ and Wa:n B ,n^ as the pre-exponents. 



Equations (13) generalize in some sense the famous bi-linear Hirota equation [Q. It can 
be obtained from (|HD in the special limit when all parameters K a:n/ 3,n 7 tends to zero. Indeed, 



let the parameters r 2jn , S3 jn and s 1 ] 1 which contains k-s tend to infinity 

l-i -2 1 

2 S ' S3 ' n ~ £ "2 



( 16 ) r 2,n ~ £ 2 +o £ 1 > S 3,n ~ S 2 ~ -£ 1 , ~ E 1 
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when e — > 0. Equating in (|H]) with a = 2, /3 = 3 and 7 = 1 the coefficients at e~ 2 and e^ 1 
we obtain respectively 

(17) T 3 n = T2 in+ e 3 , T l,n = r 2,n+ei • 

Note that in the limit flIS|) «x,n ~ K 2,n ~ ^3,n ~ £• Imposing further one more relation 
r i,n = S2,n r 3,n> one can see that the rest two equations of (|TJ]) coincide and can be written 
as a single difference equation for the Legendre variable T2 jn 

7T,nT2,n+e 1 +e 2 +e 3 T2,n = T 2i n+e 1 T2 )n+e2 +e 3 + S2, n T2,n+e 3 T2,n+e 1 +e 2 • 

In the homogeneous limit for the parameters r l n = r\ and S2,n — ^2 the latter equation 
can be rewritten in the canonical Hirota's form after obvious re-enumeration of discrete 
variables. Due to this analogy we will call sometimes the Legendre variables Tj jn as the 
triplet of r-functions. 

We conclude this section with a description of the discrete gauge invariance of the linear 
systems (|3]) and Let us require that these systems are invariant under simultaneous shift 
of the linear variables 

(18) $n !-> £ n $n , 

where the set of £ n G C. The invariance requires the corresponding change of dynamical 
variables 

£n+e 2 £n+e 2 , 

«l,n l-> 7 «l,n , ™l,n ^ ~Z ™l,n , etc. 

Sn+e 2 +e 3 Sn 

Now we have to request that gauge transformation of the parameters K a>n , a = 1, 2, 3 in 
both linear systems (^) and (f|), coincide: 

(19) « lin - f n+ ; 2+e3 Kl , n = ^ n+eien+ ; i+e2+e ^ i, n , 

Sn+e 2 Sn+e 3 Sn+ei+e 2 Sn+ei+e 3 

and similarly for the parameters K 2j n and /«3 )n . Any of this requirements leads the single 
restriction 

(20) ^n^n+ei+e 2 'Cn+ei+e 3 l Cn+e 2 +e 3 £n+ei£n+e 2 £n+e 3 'Cn+ei+e 2 +e 3 • 

The last equation can be easily solved by the ansatz 

(21) £n = ^l:n 2 ,n 3 ^2:ni,n 3 ^3:ni,n 2 



and explains appearing of all ^-dependent factors in the statement of the following Proposi- 
tion |I[ 
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4. General solution of the classical equations of motion 



As we have seen, equations ([14] ) generalize the famous Hirota equation. Actually the 
structure of ([L4]) is the same as the structure of the Hirota equation. In this section we present 
a general solution of these classical equation of motion when each three-linear equation ( 14]) 
is reduced to a pair of bi-linear equations. This form is suitable for imposing the periodic 
boundary condition. The bi-linear relations become a celebrated Fay's identities for the 0- 
functions associated with an algebraic curve of a finite genus. Periodic boundary condition 
will fix the generic algebraic curve in an unique way. The curve may depend on the type 
of the boundary conditions. We claim that for the periodic boundary conditions the found 
solution of the three-linear equations (|l4j) is the most general. We start with formulation of 
the necessary algebraic geometry objects (see e.g. f|, |5|). 

Let T be an arbitrary algebraic curve of genus g, and u = (ui,u)2, ■■■U g ) be the vector of 
g holomorphic differentials, normalized as usual: 

(22) (b io k = S kJ , (b u k = Cl 



k,j i 

J a 3 Jbj 

where Oj, bj, j = 1, . . . , g are the sets of canonical cycles on T. 
Let I: T® 2 i — > Jac(r) be the Jacobi transform 

(23) X,Y G T ^I{X,Y) = [ to GJac(r). 

Let further e (v), v G C g , e = (61,62), e% G C 9 be the thet a- function with characteristic e 
on the Jacobian Jac(r), 

(24) e (v) = ^ exp(ivr(m + ei,fim + ei) + 2m(m + ei,v + e 2 )), 

and E(X, Y) = —E(Y, X) be the prime form 1,7 6T, so that the cross-ratio 

(25) E(X,Y)E(X',Y') = £odd (I(X,F))0 £odd (I(^,r)) 
1 ' E(X,Y')E(X\Y) e £odd (I(X,Y>))e eodd (L(X',Y)) 

is well defined quasi-periodical function on T® 4 . e 0( jd is a non-singular odd theta characteristic 
such that eodd (O) = 0. We denote by 0(v) the theta-function with zero characteristic. 

There is a famous identity on T® 4 ® Jac(r), so called bi-linear Fay's identity, which can 
be written in the form: 

0(v) e(v + I(B + D,A + C)) = 

(26) =|&§8^™>°<*+^H 

where one should understand I(B + D,A + C) = l(B, A) + l(D, C) = I(A A) + I(B, C). 

Let v G Jac(r), X ni , X' ni , Y m , Y^ 2 , Z n3 , Z' n3 (n a G Z) and P, Q be arbitrary distinct 
points on the algebraic curve T. We have 
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Proposition 1. Any solution of local equations of motion ( |7^ is a particular case of the 
following: 

n,n = ^l:n 2 ,n 3 ©(In + I(<5,^m)) , 

(27) r 2 , n = & nun:i Q(I n + I(Q,Y n2 )) , 

7"3,n = s3:ni,ri2 

e(i n + i(Q,z n3 )) , 



ni— 1 U2 — 1 "3 — 1 



w/jere 

(28) in = v + K.^) + E iQ^.nj + E i(^ 3 , ^ 3 ) 

mi=0 rra2=0 m 3 =0 

an<i i/ie parameters £ and points X ni , Z' n3 enter the parameterizations of K ajn and pre- 
exponents as follows: 

£l:n 2 ,n 3 £l: n 2 +l,ra 3 +l E(Y^ 2 , Z ns ) E(Y n2 , Z' nz 



«1: 



™2,™ 3 



(29) K 2:ni 



"2 



and 



(o) 



Si: n 2 +l,n 3 sT: n 2 ,n 3 +l 


E(Y> 2 , 


Z 'n 3 ) 


E(Y n2 ,Zj ' 


^2: ni+l,n 3 ^2: ni,n 3 +l 


E(X m 


' Z n 3 


) E{X' ni ,Z' n3 ) 


S2:ni,n 3 s2:ni+l,n 3 +l 






)E(X ni ,Z'J 


S3:ni,Ti2 s3:ni+l,ri2+l 


E{X' ni 


Y 


)E(X ni ,Y> 2 ) 


S3:ni+l,n.2 S3:ni,ri2+1 


E(X ni ,Y n2 


)E{Xi nv Y> 2 ) 


si: n2,n 3 ^3: rti,ri2+l 


E{Y> 2 , 


Z n 3 ) 


E(Y n2 ,X ni ) 



Wi°i 2jn3 £,l:n 2 +l,n 3 €3:m,n 2 E(Y^ 2 , X ni ) E(Y n2 , Z n ^ 



u 2:ni,ns £l: n 2 +l,n 3 +l ^2: n\,m -^(^n 2 ' ^n 3 ) E(X ni , Z£ 



(30) 



"3 ' 



"l:n 2 ,n 3 



W 3:ni,n 2 ^2: m+l,n 3 4^3: ni,n 2 +l E (Y^ 2 , X' ni ) E (X ni , Z, 



4 r i:n2+l,n 3 6:ni,n 3 +l EiY^ Z' n J E(X ni , Z n J 



n 3 i 



w£> ni m 6:711,713 6:ni+l,n 2 +l ^(^4 ' ^m) ^ (-^ni ' ^nj 



Proof: In order to prove that substitution (|IB|) with r-functions given by (ETO solves ( II ) 
one should use repeatedly the Fay identity. For example, equation ( |l4l) for the choice at = 1, 
(5 — 2 and 7 = 3 is the consequence of two Fay identities (p6|) taken for two sets of divisors 
(A 5, C, D) and (A', B\ C, D') with identification A — A' — X ni , B — B' — Y n2 , D — D' — 
Y^ 2 , C = Z' nz and C = Z n3 . The appearing of the ratios containing parameters £ a :np,n-, 
are due to gauge invariance of the linear system (§) or @. Gauge parameters s r a:n /3 ,n 7 , the 
divisors X ni , X' etc. as well as period matrix Q and the point on the Jacobian v are free 
parameters of the general solution. By imposing the periodic boundary condition (which 
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means in particular the fixing of the size of the system) they may be related to the values 
of the integrals of motion. ■ 

One may solve expressions (|30D in order to avoid the rations of the pre-exponents. To do 
this we need to introduce several extra parameters A ni , B n2 , C m . Then, one obtains for the 
"first" block: 



£l:n 2 +l,n 3 E[Y. ' , Z ns ) E(C n3 , Z n J T 2>n 



(31) w hl 



Ci:n 2 +i,n 3 +i E{Y^ 2 ,Z n3 ) E(C n3 , Z' n3 ) T 2in+e3 

6:n 2 +l,n 3 E ( Z n 3 > Y n 2 ) E ( B n 2 ' Y n 2 ) T3,n+e 2 
£l:n 2 ,n 3 E \ Z m^nJ E i B n 2 -> Y n 2 ) r 3,n 

_ _ £l:n 2 ,n 3 £l:n 2 +l,n 3 +l E ( Y n 2 ' Z n 3 ) E ^Xn 2 i Z 'n z i 
£l:n 2 +l,n 3 £l:n 2 ,n 3 +l E ( Y n 2 ) Z 'n 3 ) E \¥n 2 > Z n 3 , 



for the "second" block: 



u ^ _ t,2:m,n 3 E{X nv Z' m ) E(C n3 , Z n3 ) Ti ; n 

E{X ni ,Z n3 ) E(C n3 ,Z^ 3 ) Ti jn -\ 



he :1 



w — _ ^2:m,ra 3 E(Z n:j ,X ni ) E(A ni , X ni ) r 3>n 

6:ni+l,n 3 E(Z n3 ,X ni ) ^(Au'^ni) r 3,n+ei 



^ _ _ ^2:m+l,n 3 ^2:m,n 3 +l E(X ni , Z n J E(X' ni , Z' n3/ 

^2:m,n 3 ^2:ni+l,n 3 +l E{X' nv Z n ^) E{X n ^ Z' n J 



and the "third" block finally: 

«3,n 



^3:ni,n 2 +l E(X^ Y n J E(B n2 ,Y^ 2 ) Ti >n+e2 
6:ni,n 2 ^(^nii^W E(B n2 ,Y n2 ) Ti, n 



(33) u> 3 ,n 



^3:m+l,n 2 +l ^O^j^ni) ^(Au ' ^nj r 2,n+ei 



_ ^3:ni,n 2 ^3: ni+l,n 2 +l E(X' ni , Y^J E(X ni , Y^ 
C3:ni+l,n 2 ^3:ni,n 2 +l E(X ni , Y n2 ) E(X' ni ,Y^ 



where r ajn are given by (|27j ). In addition, the discrete Baker- Akhiezer function <3> n , obeying 
the whole set of (§) (and therefore, (f|)) and normalized by $o = £o, is given by 



(34) 
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where 



(35) 



*(o) ( p \ TT E ( P ^ X m 1 )E(A mi ,X' mi ) 

n " li E{P,X> mi )E{A m ^X mi ) 

"fy E(P,Y m2 )E(B m2 ,Y^ 2 ) "j-j- 1 E(P, Z ma )E(C m3 , Z' m3 , 

X 



R(P V \E( R V \ 11 



:Q E(P, Y; i2 )E(B m2 , Y m2 ) A £(P, Z; i3 )E(L7 m3 , ^ 3 ) 

5. Finite lattice with open boundary 

As one could see, all the considerations until now were local and did not take into account 
the size and the boundary of the cubic lattice. To describe the global characteristic, such as 
integrals of motion, we should specify this external data. 

The most natural is the case of the finite cubic lattice with open boundary conditions and 
the corresponding Cauchy problem. For the cubic lattice of the size Ni x N2 x N 3 we fix the 
coordinates ni,ri2,n 3 of 3d vector n (|l|) as follows: 

(36) < n a < N a , a = 1,2, 3. 

A = N1N2 + N2N3 + N3N1 edges incoming to the cubic lattice correspond to the 2A initial 
data 

^1:712, ri3 = ^I,0ei+n2e2+ri3e3 > "^1:712,713 = "^l,0ei +71262+71363 j 

(37) 

^2:77,1,713 — ^2,niei+0e2+n 3 e3 1 W2: ni ,n 3 — W2,niei+0e2+77,3e 3 5 
:ni,n2 

= U 3 

,niei+n2e2+0e3 > "^3:711, J12 — 

W 3 

,niei+ri2e2+0e3 • 

Equations of motion @-(||), being applied recursively, defines the transformation from the 
initial data fl3"T|) to the 2A final data 

M l:ri2,ra3 = ^l,iViei+n2e2+n3e3 > W l:n2,nz = W l,N 1 ei+n2e2+me 3 1 
(38) U 2:ni,ri3 = u 2,niei+N2e2+n3es 1 ^2:n\,nz = ^2,niei+A'2e2+n3e3 j 

U 3:n 1 ,n2 = U 3,niei+n2e2+Nze3 1 w 3 ;n\,n2 = W 3,n 1 ei+n2e2+Nse3 ■ 



Suppose that the initial data fl37|) are generic. The natural question arises. How one can 
invert the parameterizations described by the formulas (|31j)-(|33|) in order to restore the 
algebraic geometry data in terms of (|37"D and the parameters K a:n/3 ,n 7 (the total number of 
parameters is equal to 3A). 

In general the solution of this problem is not unique. Evidently, the same initial data 
may be parameterized by the infinite number of the different sets of algebraic geometry data 
with sufficiently high genus. But there exists one preferred compact Riemann surface, which 
is minimal in the set of all possible curves parameterizing the dynamics (^-§). 

In the Cauchy problem with (|37j) and (j3~8|) this curve appears as follows. Let the initial 
data are parameterized in the terms of some F One may write 
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Then due to (^7],^) the solution of the Cauchy problem is 

M a:n (3 ,n 7 = M o:n (3 ,n 7 ( v + T Q ) , W a]n/3jni " , o:n ( j,n 1 ( v "I" T Q ) , 

where 

ATi N 2 N 3 

(39) T x = ^I(X^,X ni ), T 2 = X) 1(1^, yj, T 3 = 53l(Z^,Zj. 

ni=0 ri2=0 n,3=0 

The transformation from the initial data to the final ones is the evolution if 

(40) Ti = T 2 = T 3 = T mod Z 9 + VLZ 9 . 

In what follows, all the relations of the type ( (4"0"| ) will be understood modulo Z 9 + QZ g . 

The evolution conditions ( fTO|) mean that due to the Abel theorem there exist two mero- 
morphic functions A(P) and /i(-P), P G T, with the divisors 

(-^) = 53 ~ 53 ~ 53 + 53 ' 

, nieZ^ nieZjv-,^ n3GZjv 3 n3GZjv 3 

(a*) = E F - - 53 K- 53 z «.+ E ^ 3 • 

ri2&N 2 "2GZjv 2 "3SZjv g W3£ZjV 3 

This means that T is the compact Riemann surface given by a polynomial equation 
Ja(A, /i) = (see e.g. Theorem 10-23 in |J). Moreover, (|4l|) fix the generic structure 



of Ja(A, /i) (later we will give it, see (|4q) ). The key observation is that due to the evolu- 
tion conditions (|40| ) Ja(A, //) as a functional of the initial data should produce the set of 
invariants of the evolution, so the curve is the spectral one. 

Ja may be derived with the help of the linear system of the type (|j) written for the whole 
auxiliary plane. Let us fix the position of the auxiliary planes corresponding to the initial and 
final data. The auxiliary plane for the initial data crosses all the incoming edges of the cubic 
lattice while for the final data it intersects all the outgoing edges. These auxiliary planes 
plays the role of the two-dimensional space-like surfaces and discrete evolution corresponds to 
the translation of the auxiliary plane into direction perpendicular to this space-like surfaces. 

All the objects, namely the dynamical and auxiliary linear variables on the space-like sur- 
face, can be numbered by the two-dimensional discrete index. Let us choose the numbering 
for the linear variables in the form similar to numbering of the initial (|37D or final (|38D data: 



< ^ > 0ei+n 2 e 2 +n3e3 ( ^ > l:n, 2 ,n 3 j ^ > niei+0e 2 +n 3 e3 < § > 2:ni,n 3 j 

(42) 

^mei+n2e2+0e3 ^3:ni,ri2 ? ^ — — i ) ^ 1, 2, 3 . 



An example of such enumeration on the auxiliary plane in the simple case N\ = N 2 = N 3 = 2 
is shown in Fig. [5]. 
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m=0 



ni=l 




2:22 



n 3 =0 n 3 =l 

FIGURE 5. Example of enumeration of the linear 



■ variables. 



Let C(X, ji) be the matrix of the linear system 



(43) 



— $i : „ 2 +i in 3 — <&l- n2+ l yn3+ \Ul 

:n 2 ,n 3 :n 2 ,n 3 

^2:ni, n 3 ( ^ ) 2:?ii, n3+1^2:ni, «3 ^2:ni+l,n 3 ^2:ni;n 3 ~^~ 

+ &2:ni+l,n 3 +l K '2:n 1 ,n 3 u 2:n 1 ,n 3 w 2:n 1 ,n 3 j 
= $3:ni,ri 2 +l — < ^3:n 1 ,n 2 M 3:rii,n 2 + &3:n 1 +l,n 2 +l w 3:n 1 ,n 2 + 
+ ^3:ni+l,n 2 K '3:n 1 ,n2 U 3:n 1 ,n2 W 3:n 1 ,n2 

written in the matrix form = $-£(A, fi), where the linear variables satisfy the identification 
conditions 

(44) $ 1:0 ,n 3 = $2:0,n 3 , $1:712,0 = $3:0,n 2 , $2:ni,0 = $3:ni,0 , 

and the quasi-periodicity conditions 

$3:^,712 _ ^1:^,^3 ®3:ni,N 2 _ $2:ni,N 3 ®l:N 2 ,n 3 _ ^2:JVi,n 3 



(45) 



x z y z y 

for the boundary domains on the auxiliary plane. Parameters 

y 



X 



(46) 



A 



x 



z 



[I 



z 
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FIGURE 6. Newton's polygon. 



are the complex numbers and we call them the spectral parameters. Due to ( f44|) and 
it is clear that the total number of the independent linear variables in the system (|43D is A 
and C(X, fi) is A x A square matrix. 
Define 

(47) J A (A,/x) = det£(A,/x) I ]J[ u 1:n2>r 

\ri2,m 

(fi?]) is a normalized Laurent polynomial {Jo,o = 1) of the spectral parameters A and fi, 

{0 < a < N 2 + N 3 
< b < N, + N 3 . 
-Ni <a-b<N 2 

Domain II (Newton's polygon of Ja) is shown in Fig. |6]. One may prove (see jl]) that the 
coefficients of this Laurent polynomial are invariants of the evolution, i.e they are the same 
for the initial fl37f ) and final fl38f ) data and can be served as a source of the algebraic geometry 
data for the solution of these equations with open boundary conditions. 

The requirement that the linear system (^) has non-trivial solution is equivalent to the 
fact that the spectral parameters belongs to the algebraic curve: 

(49) F=(A,//)er A ^ J A (X,fi) = 0. 

Assuming that all the incoming data together with parameters K a:n/3 ^ are in general position, 
one may calculate the genus of the curve ( ^9|) using Newton's polygon. For the cubic lattice 



of the size N± x N 2 x A^3 such that N± > N 2 > A^3 Newton's polygon associated with this 
curve (f49D is shown in Fig. || The genus of the curve Ta is equal to the number of internal 
points of this polygon 

(50) g A = (iViiV 2 + N1N3 + N 2 N 3 ) - (A^ + N 2 + N 3 ) + 1 . 

g coefficients J a ^ in the Laurent polynomial J/\(A, /i) corresponding to the internal points of 
the Newton's polygon are related to the moduli of the algebraic curve and so give the 
period matrix Coefficients J a ,b corresponding to the perimeter of the polygon are related 
to the divisors of the meromorphic on Ta functions A,/i and so give the set of X , Z' n3 . 
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From (|HD and (^) we may get explicitly 



N fr' E(P,X ni )E(A ni ,X' ni ) E(P, Y; 2 )E(B n2 , Y n 
11 R(P Y< \TP.( A.. Y \ 11 



(51) 



; o E(P,X^)E(A ni ,X ni ) 11 E(P } Y n2 )E(B n21 Y> 2 ) 

^ E(P,Y n2 )E(B n2 ,Y> 2 ) ^ J5(P, Z'JE{C W Z n3 ) 
11 R/PVME/R V wl 



£(P, Y> 2 )E(B n2 , Y n2 ) 11 £(P, Z n3 )E(C7 n3 , 



Previously in [[Jj] a bit different approach was used. Due to the Theorem 2 of [0], the 
solution $(A, //) of the linear problem $ • £(\, fi) = as the vector of meromorphic functions 
on the curve Ta is given by (|34|) on Jac(rA). With the help of (|3^ ) the parameterizations 
( |27| - |30D may be restored uniquely 

We conclude the discussion on the Cauchy problem by the calculation of the degrees of 
freedom of the algebraic geometry parameterizations. With A = N1N2+ N2M3+N3N1, A' = 
N1 + N2 + N 3 and g — A — A' + 1, we have g moduli of Ta, g complex numbers v 6 Jac(rA), 
g independent cross-ratios of £ n , 2 A' — 2 independent divisors X ni , Z' n3 and A' arbitrary 
divisors A ni , B n2 , C n3 minus one degree of freedom corresponding to the arbitrariness of Q: 
totally 3A parameters. It proves the completeness of the parameterizations in the terms of 

r A . 



6. Backlund transformation 

Curve Ta corresponds to the generic position of the initial data and open boundary con- 
ditions. For the application to the integrable spin models the periodical boundary conditions 
T = 0, i.e. 

(52) ^K^XJ = ^K^nJ = E ^ns^ns) = 

niSZjVj «2GZjv 2 nz&M 3 

are important. Periodical boundary conditions reduces the spectral curve Namely, 
due to fl52p, the parameters x,y,z in ( [I5D become the meromorphic functions with divisors 
(x) = J2 ni X ni — J2 ni X' ni an d O). Since now x, y, z are meromorphic functions (not only 
their ratios) the curve can be defined by an algebraic equations for any pair from x, y, z: 

(53) Ji(y,z) = J 2 {x,z) = J 3 (x,y) = 0. 

The Laurent polynomial J&(x/z,y/z) is not a generic one and Ja = is a consequence of 
the algebraic relations (|53|). 

In this section we describe explicitly the reduced spectral curve and the meromorphic 
functions, which uniformize it for another choice of the auxiliary plane. The interpretation 
of the problem now differs from the Cauchy one for the open cubic lattice. The discrete 
evolution can be identified now with a sequence of the Backlund transforms for the square 
auxiliary plane. 
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First we define the new position of the auxiliary plane. Let it crosses the ingoing edges 
for the vertices of the 3-dimensional lattice with coordinates n2&2 + ^363. With this choice 
the role of the variables u a and w a , a = 1,2,3 becomes different: U2, a , ii>2,n, U3,n, u>3 in 
are auxiliary, while Ui }Il , Wi >n are dynamical. The evolution will describe the change of these 
dynamical variables with respect to the "time" n\ with periodic boundary conditions implied 
in the second and third directions: 

(54) Mi 

:n 2 ,n 3 +N 3 = Ui 

:n 2 ,n 3 , :n2+N2,ri3 ~ W l:n 2 ,n 3 +N 3 — Wl:n 2 ,n 3 . 

Of course, the auxiliary for this evolution variables are also different from layer to layer 
according to the discrete equations (0)-®. We will keep the open boundary conditions in 
the first direction. 

Consider one layer of the cubic lattice corresponding to n\ = 0. Equations of motion @ 
and the periodical boundary condition (|5i]) yields an implicit transformation 

(55) Ul:n 2 ,n 3 ,Wl:n 2 ,n 3 > u l:n 2 ,n 3 5 W l:n 2 ,n 3 ■ 

According to (p^-|30|), one may define 



(56) ^I:ri2,n3 — ^l:n 2 ,n 3 ( v ) , ^l:n 2 ,n 3 

where v e Jac(r) for some algebraic curve T. Then the transformation fl5"5"|) can be written 
in the form 

(57) <n 2 ,n 3 = M l:n 2 ,n 3 (v + I(X' Q , X )) , 11]'^ ^ = W lin2>n3 (v + I(X' Q , X )) . 

The periodicity conditions (|54]) in the parameterizations ( j27| - |3(i| ) means J2n 2 ez I(^n 2 ^n 2 ) = 
Tlin 3 & N I(Z' n3 , Z n3 ) = 0, (see (|52"D). Repeating the consideration from page |12|, we may 
conclude that there exist two meromorphic functions y and z on V such that 

(58) (V) = £X ~ K , (*) = E Z n 3 ~ Z n 3 > 

n 2 n 3 

and a compact Riemann surface T may be defined by a polynomial equation Ju{Vi z) — 0. 
Moreover, the structure of (y) and (z) fixes the algebraic form of F (see (|6i~|) later). 

We call transformation ( |55|) the Backlund transformation for the two-dimensional square 
lattice, since fl55|) is a canonical transformation and it conserve the integrals of motion. In 
order to specify the integrals of motion we again consider the linear system for the chosen 
auxiliary plane. This linear system is a system of N2N3 equations 

(59) j n2jH3 = 0, $_ 1>n3 = y _1 $jv 2 -l,n3 > ®n 2 ,N 3 = Z$n 2 ,0 , $-l,JV 3 = y _1 ^iV2-i,o , 

where < 112 < N2, < < N% and linear form j n2 ,n 3 is defined by the first line in (|3]) with 
slightly shifted enumeration of the auxiliary linear variables $ n2 ,n 3 

^gQ^ jn 2 ,n 3 = ^n 2 ,n 3 ^n 2 ,n 3 +\U\: n 2 ,n 3 + ^n 2 — l,n 3 W\; n 2 ,n 3 + 

^n 2 — ljTia+l^li n 2 ,n 3 'U'l: n 2 ,n 3 ^l: n 2 ,n 3 • 
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Equations in ( |5PD which includes the spectral parameters y and z describe quasi-periodical 
boundary conditions for the linear variables $n 2 ,n 3 - 

Let L(y, z) be the complete matrix of the coefficients of the system fl59f) , j = $ • L(y, z). 
Define a Laurent polynomial Jn(y, z ) 

N 3 N 2 

(61) J n (y, z) = det L(y, z) = J2/2 J ^ b V* z " ' J °>° = 1 ' 

a=0 6=0 

J a ^ in this case are the invariants of the Backlund transform (|55|), see HJ. In order the 
system $(P) • L(y, z) = has a nontrivial solution, the spectral parameters y and z should 
belong to the spectral curve 

(62) P=(y,z)eT n & J D (y,z) = 0. 

Assuming again that the dynamical variables U\ :n ^ nz , W\; n2jn3 and parameters ni :n2 ,n 3 are 
generic, the genus of the spectral curve is equal to g\j = (N 2 — 1)(N 3 — 1). Recall, 
the form ( |BTD of the polynomial Jq follows uniquely from the conditions (y) = Y n2 — Y^ 2 , 

In the backward direction, Theorem 2 from says that the non-normalized solution of 
$(P) ■ L(y, z) = as the meromorphic function of T\j is given by 



E(P, Y m2 ) E(B m2 , Y^ 2 ) -r-j E(P, Z m3 ) E {C mz , Z' mz ) 



^n 2 ,nA r ) 11 prp V , \ rp( id V \ 11 



X 



„,, :[] E(pY^)E(B m2 ,Y m2 ) ^i E(pZ^)E(C m3 ,Z m3 [ 

V ^) / n 2 -l n 3 -l 

xe v + i(q,p)+^ i(y; 2 ,y m2 ) + ]T i(z m3 ,z m; 

\ m 2 =0 m3=0 

where theta-functions are constructed by means of the period matrix fin of the algebraic 
curve rrj. Using formula (|63| ) we may find the corresponding expressions for the dynamical 
variables ui :n2j1l3 and wi- n2tns to observe that they coincide with those given by the Propo- 
sition |I] and arbitrary curve T being identified with the spectral curve Explicit form of 
the spectral parameters uniformizing the spectral curve ]?□ are 



IT E(P,Z n3 ) EjC^KJ 



(64) 



v(p) TT E(P,Y n2 ) E(B n2 ,Y; 2 ) 

JJ E(P,Y> 2 ) E(B n2 ,Y n2 ) 



In case when the algebraic curve used for the construction of the general solution becomes 
a rational one, the one step of this evolution n\ i— > n\ + 1 with open b.c. in the 1-st direction 
can be identified with creation of additional soliton. This will be explicitly demonstrated in 
the next section. 

In this Section we have established the origin of the curve Tu appearing when the peri- 
odical b.c. are imposed in 2-nd and 3-rd directions. Evidently, Jn(y, z) should be identified 
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with Ji(y, z) in fl5"3|) . In the same way the periodical b.c. may be imposed in any other pair 
of directions, and the corresponding J 2 and J3 also are the spectral determinants. If the 
periodical b.c. are imposed in all three directions, then any of three relations J% = 0, J 2 = 
and J3 = should define the same curve, and since the principle of generic data in this case 
has been lost, the genus of the curve may be established as 

(65) g < min of (N, - 1)(N 2 - 1) , (N 2 - 1)(N 3 - 1) , (JV 3 - l)(JV a - 1) . 



7. Rational limit 

Let us again forget on the boundary conditions and consider a rational limit of the alge- 
braic geometry solutions to the discrete equation of motion This rational limit for the 
0-function (121) on a Jacobian of an algebraic curve corresponds to (see |H) 



and 

(67) » 0, e-^" 1 — > ^-^) (Pk-Pn) dg 4 ^ 

(?fe - Vn) \JPk ~ Qn) 

Thus the set of c4, n is the reminder of the period matrix, while the set of f n is the reminder 
of v. The prime forms in the rational limit are the prime forms on the sphere: 

E(A,B)E(C,D) (A-B)(C-D) 
{ ' E(A,D)E(C,B) (A-D)(C-B)' 

For g — 1, 2, . . . and the sequence of the parameters 

(69) P0,g0,/0; ■■■ Pg-l,qg-l,fg-l = {Pk, <?fc, /fc}fc=J 

we introduce the rational limit of the O-function (see the appendix of ||) 

(70) H(°\{p k , Qk , f k }lJ ) = . 

life - ?i) 

i>j 

We set = 1 by definition. Note that if all parameters /& vanish H^({p k , q k , 0}^~q) = 1 
as well. Let the function a k (z) be 

(71) a„{z) = ^— - • 

qk - z 

When g = the equations (0) in the rational limit have the simple solution r a n = 1 due 
to identity 

(72) r a = 1 + S/3 + s' 1 , 
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where 

Xa X X Xfj X a X'a 



Xp — X 7 X a — X 7 X a — x^ 

(73) 

X a — X 7 X a — Xp 



X a — Xp X' a — X 7 

and (a, /3, 7) is any even permutation of the set (1,2,3). The notations in (ff2j) and (f73|) are 
related to those in (|ToD and the Proposition [I] as follows: r a = r ajn , s a = s a>Il ] X\ = X ni , 
X[ = X' ni , X 2 = Y n2 , X' 2 = Y; 2 , X 3 = Z nz and finally X> = Z' nr 



Proposition 2. The soliton solutions of the equation (~T$) are given by the formulas 



Tl,„ 



H (9) 



In: 



n: k 



o~k{X ni ) J fc=0 



9-1 



(74) r 2 , n = 



n: fc 



o- k (Y r 



9-1 



T 3,i 

where I n - k , ni,n 2 ,n 3 G Z 





r 


[(T k (Z n3 ) t 


fe=0/ 



i=0 a k( X mi)J Vm 2 =0 ^vj) Vrnai 

and parameters r Q n and s ain are defined by the formulas (|75|j. 

The proof of this Proposition is based on the rational variant of the Fay's identity 

(A - D)(C - B) H^({f k ^§-}U H^\{f k a -^hu) + 

(76) (A -B)(D- C) H^\{f k ^§}i =1 ) H^\{f k ^Y k=1 ) = 

(A -C){D- B) H^({f k }U H^\{f k ak ^\ a f \ y k=1 ) 

a k {B)a k {D) 

described in ||. 

So far the parameters p k , q k and f k (|69| ) are arbitrary complex parameters and the solution 
given by the Proposition is relevant for lattice infinite in all directions. Now we would like 



to impose periodic boundary conditions (53) in directions 2 and 3 and interpret the evolution 



along the direction 1 as a sort of Backlund transformation which create the solitons. To do 
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this we first note that boundary conditions ([3]) are equivalent to the following algebraic 
relations for the parameters p and q: 

(77) fr 1 ^ = rr 1 ^) = i. 

rt2=0 113=0 v °' 

One may verify that for the parameters Y n2 , . . . , Z' n3 in general position the system of equa- 
tions (|77|) has exactly g = (N 2 — 1)(N 3 — 1) non-equivalent solutions (equivalence means 
that if (p,q) is a solution of (|77|) , then (g,p) is also the solution). Let us choose this set of 
solutions as sequence ( |69| ) leaving parameters f k , k = 1, . . . , g to be free. 
Using this freedom let us redefine the 'amplitudes' f k as follows: 

(78) fk = F k -a k (X k ), 

where we recall that parameters p k and q k of the functions ([fl]) are already fixed by the 
system of equations ([77]), while parameters X k are still free. Let us consider the solutions 
for the tau-functions r 2 n and r 3 n given by the Proposition |2| and with redefined amplitudes 
([78|) at the value of the discrete coordinate n\ — 0. Using the freedom in parameters X k we 
send 

(79) X k ^p k . 

It is clear that in this limit all components of the set {J n: k } 9 k =o vanish and according to 
definition ([70]) we will have r 2jn = r 3 n = 1. In other words, we obtained for these tau- 
functions at ri\ = a zero-soliton solution. Let us repeat the same procedure at ri\ = 1. 
Namely, again consider generic solution for these tau-functions given by the Proposition ^| 
and then make a limit ([791). It is clear that now in the set {I n :k} k Zo one element which 
correspond to k — will survive and as result we obtain one-soliton solution. Increasing 
ni results in increasing the number of solitons. One may see that the maximal number of 
solitons which can be reached by this procedure is equal to g = (N 2 — 1)(N 3 — 1). See detailed 
description of this phenomena in the simplest situation discussed in the paper ||. 

But before conclude this section we would like to give one more explanation of the solitons 
creation procedure. Let us consider the equation (|H]) at rii = for a, (3, 7 = 1,2,3 and 



homogeneous or zero-soliton r-functions T2°n 2 ,n 3 = T 3?n 2 ,n 3 = 1- This is a difference with 
respect to "space" coordinates n 2 and n 3 linear equation for the function ri :n2 n3 . Using 
simple algebra one may verify that besides trivial solution T 1:n2 n3 = 1 this equation has 
a solution of the form T\- nnn „ = ni 2_ n — tt^t ni 3_1 n — 7^4, where parameter p of the 
function ctq is identified with X . The complete solutions of linear difference equations is 



1 ff o(Q "rj 1 a v{ z 'm 

m,2=0 m3=0 ° 



with arbitrary F , Now solving the equation (P^4|) for a, /3, 7 = 2, 3, 1 and 3, 1, 2 with found 
r i:rL n 3 we can fi nc ^ ^ ne va hies of the r-functions T2:n 2 ,n3 and T 3:n2)n3 at the discrete time 
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rii = 1. They will have similar to ( |50"D one-soliton form. Using these solutions again in the 
equation ( fH]) with a, (3, 7 = 1,2,3 we will find two-soliton solution for the function 7i :n2)n3 
and then two-soliton solutions for the functions T 2:n2in3 and T 3:n2in3 at the next value of the 
discrete time n 2 = 2. It is clear that this procedure can be continued and demonstrate the 
equivalence of the discrete evolution parameter rii to the number of solitons. This simple 
explanation finally justify that investigated in this paper the discrete dynamic given by the 
equations of motion is a set of consecutive Biicklund transformations. 

8. Discussion 

As it is well known [T|], the dynamics of parameters of 3d spin models is equivalent to the 
dynamics (H-§). So the classical solutions discussed in this paper can be used for the spin 
models of different types. The solutions with periodic boundary conditions can be utilized 
for construction generic spin model, such that their Boltzmann weights are parameterized 
by the theta-functions of the higher genus. The integrability of such spin models is based 
on the modified tetrahedron equation ||3_ 0f| . This should include the different generalization 



of the famous Chiral Potts model. 

On the other hand the solitonic solutions are convenient for completely inhomogeneous 
ZBB model and open a way to develop the quantum separation of variables for ZBB model || . 
In particular, these solutions with parameters fk 7^ allows one to construct the complete 
family of isospectral deformations of ZBB transfer-matrix (see also || for realization of this 
program in case of relativistic Toda chain model with spin degrees of freedom). 
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